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ABSTRACT 

This  report  constitutes  Part  2 of  a two  part 
report.  Part  1 "BEAR  Buoy:  The  Engineering  Documenta- 
tion for  Scientific  Application",  dealt  with  details 
of  the  system  design  and  calibration.  This  report 
first  develops  the  mathematical  model  used  to  determine 
the  mooring  configuration  for  varying  conditions  of 
external  forcing.  Second,  the  method  used  to  correct 
sensors  (thermistors)  for  vertical  displacements  is 
I presented.  Included  are  case  study  comparisons  of  the 

I model  predictions  vs  actual  measurements. 
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1 . 0 Introduction  and  Background 

Oceanographic  measurements  from  moored  buoys  are  con- 
taminated by  mooring  motion.  Previous  studies  (eg,  Wunsch 
and  Dahlen,  1974)  have  shown  that  deep  instrumentation 
located  on  surface  buoyed  mooring  lines  can  undergo  vertical 
excursions  as  large  as  several  hundred  meters.  Large  changes 
in  sensor  depth  result  from  low  frequency  horizontal  meander- 
ings  of  the  surface  buoy  and  changes  in  configuration  of 
the  mooring  line  under  the  action  of  the  surface  forces  of 
wind,  wave,  and  current  and  subsurface  drag  due  to  current. 
High  frequency  motions  are  present  as  well  due  to  the  motion 
of  the  surface  wave  field  and  due  to  line  strum  on  the 
mooring  line  caused  by  vortex  shedding.  For  an  operational 
moored  system  a model  is  needed  which  can  predict  the  motion 
and  configuration  of  the  mooring  under  varying  conditions 
of  external  forcing.  Given  the  predicted  configuration 
the  measurements  can  then  be  corrected  for  vertical  dis- 
placements . 

This  report  constitutes  Part  2 of  a two  part  report. 

Part  1 "BEAR  Buoy;  The  Engineering  Documentation  for 
Scientific  Application",  dealt  with  the  details  of  the 
system  design  and  calibration.  Part  2 details  the  ana- 
lytical model  used  to  predict  mooring  configuration  and 
the  auxilliary  model  used  to  correct  the  sensor  signal 
for  vertical  displacements.  Included  are  case  study 
comparisons  of  the  model  predictions  vs  actual  measure- 
ments. These  were  used  to  validate  the  model  and  to 
establish  the  limits  of  ambiguity  for  the  displacement 
corrections . 
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1 . 1 The  BEAR  Buoy  System 

The  buoy  was  fabricated  from  two  NAVAIR  cruiser 
mooring  buoys.  The  two  buoys  were  strengthened  and 
welded  together  to  provide  the  internal  volume  for 
instrumentation  and  additional  buoyancy  to  support  the 
thermistor  string  and  mooring.  The  buoy  has  a total 
weight  of  13061  Kg  and  is  filled  with  expanded  poly- 
urethane foam  which  provides  a positive  buoyancy  of 
10715  Kg  per  meter  of  buoy  draft. 

Buoy  motion  was  measured  and  recorded  by  a 
separate  on-board  recording  data  system  supplied  by 
the  National  Data  Buoy  Office,  Bay  St.  Louis, 
Mississippi.  The  system  - Portable  Ocean  Platform 
Motion  Instrumentation  Package  (POPMIP)-  is  battery 
operated  and  self-contained.  The  POPMIP  remains  in  a 
stand-by  mode  until  activated  by  radio  command.  After 
turn-on,  the  system  sequentially  samples  each  of  6 
analog  sensors  3.82  times  per  second  over  a turn-on 
interval  of  2 minutes  cind  then  returns  to  the  stand-by 
mode.  Each  sample  is  converted  to  a three-digit  BCD 
word  cind  recorded  on  1/4  inch  magnetic  tape.  The  in- 
ternal format  provides  a total  of  90  independent  samples 
per  tape. 

The  mooring  sensors  included  the  tensiometer, 
inclinometers,  and  pressure  gauge.  The  strain  gauge 
tensiometer  was  mounted  in  the  clevis  at  the  base  of 
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the  buoy.  The  tensiometer  provided  a measure  of  the  sum 
of  the  forces  acting  on  the  buoy  due  to  the  surface  forces 
of  wind,  current,  and  waves,  to  the  sum  of  the  subsurface 
drag,  and  to  the  weight  of  the  mooring.  Two  inclinometers 
were  mounted  on  the  thermistor  cable:  one  near  the  surface 
and  one  at  the  lower  end  of  the  thermistor  string.  The 
inclinometers  were  filled  with  silicon  oil  of  60,000 
centistoke  viscosity  in  the  upper  unit  and  30,000  centi- 
stoke  in  the  lower,  resulting  in  a time  constant  of  ap- 
proximately 60  seconds  for  both  units.  This  served  to 
damp  out  the  major  portion  of  the  higher  frequency  cable 
vibration.  A pressure  gauge  located  at  a depth  of  650 
meters  flooded  during  the  buoy  deployment  and  was  not 
operational.  The  information  from  the  strain  gauge 
tensiometer  and  the  two  inclinometers  provide  inputs 
to  an  analytical  model  which  was  used  to  estimate  the 
thermistor  depth. 

The  mooring  system  as  shown  in  Fig.  1 was  comprised 
of  the  thermistor  cable,  nylon  mooring  line,  and  the 
ground  tackle  (chain  sections  and  anchor) . Table  1 sum- 
marizes the  lengths  of  each  section  under  no-load  and 
load.  The  static  load  includes  the  stretch  of  each 
member  computed  for  conditions  of  zero  surface  and  sub- 
surface forcing. 


Table  1 


BEAR  Buoy  Mooring  Components 


f 

Section  ■ 

Length  (m) 
Unstretchrdi Static 

Load 

Percent 

Stretch 

wt/unit 
(Kg  m“ 

The  thermistor  cable  comprised  the  upper  portion  of  the 
mooring  system  and  consists  of  an  inner  core  of  37 , number  18 
conductors  in  a wet  core  configuration.  It  is  protected  by 
a double  lay  torque  balanced  steel  armor  cable  with  a rated 
breaking  strength  of  33,000  kq.  The  armored  cable  was 
mechanically  connected  to  the  buoy  through  a clevis  and 
then  led  through  a pipe  to  the  top  of  the  buoy  where  it 
was  connected  by  waterproof  connectors  (Fig.  1) . Plastic 
ribbon  fairing  was  interwoven  into  the  armor  to  a depth 
of  1000  meters  to  reduce  the  strumming  of  the  cable. 

The  lower  end  of  the  thermistor  string  was  connected 
to  a 2-inch  nylon  mooring  line  with  a rated  breaking 
strength  of  41,000  kg. 

The  lower  end  of  the  nylon  mooring  line  terminated 
in  5 shots  (27.4  m/shot)  of  1 1/8  in.  stud  link  and  14  shots 
of  l>s  in.  buoy  chain  weighing  a total  of  12020.1  kg.  The 
chain  acted  as  a varying  dynamic  load  for  the  mooring  and 
also  served  to  virtually  increase  the  scope  of  the  anchor 
to  maximize  its  holding  power.  The  anchor  was  an  1800  kg 
Navy-type,  with  an  expected  holding  power  of  5400  kg. 


BUOY  DIAMETER 
3.6  m. 


/"  THERMISTOR  CABLE 
(15  THERMISTORS  = 
T2-TI6) 


TENSIOMETER, 

INCLINOMETER, 

THERMISTOR-TI 


1935  m. 


_ II  CLINOMETER 


2“  DIAM.  NYLON 


4795  m. 


2820m. 


1800  Kg 
ANCHOR 


19  SHOTS  OF 

CHAIN 

(l•/e•l•/2) 


BEAR  BUOY  SYSTEM 


Fig.  1 
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1 . 2 Background  for  the  Development  of  the  Mooring  Model 

The  analysis  of  a single  point  mooring  can  be  divided 
into  two  problem  areas.  1)  The  static  problem  which  is 
concerned  with  forces  acting  on  the  system  - wind,  wave, 
and  current  - that  are  considered  steady  over  some  discrete 
time  interval.  This  has  been  treated  extensively  by  Pode 
(1951),  Wilson  (1961),  and  Morrow  and  Chang  (1967). 

2)  The  second  problem  is  that  of  the  dynamics  which 
treats  oscillatory  wave  forces  and  impulse  forces  in 
wind  and  current.  The  latter  problem  has  been  treated 
extensively  in  the  past.  For  example: 

a)  the  dynamic  response  of  an  inelastic  line  in 
which  the  motion  at  the  end  points  was  assumed 
(eg;  Walton  and  Polachek,  1959), 

b)  deep  water  flexible  moorings  with  the  surface 
float  (or  ship)  subject  to  harmonic  waves  (eg; 
Wilson  and  Garbaccio,  1967;  Langer,  1964),  and 

c)  three  dimensional  mooring  dynamics  in  which  the 
line  hysteresis  was  considered  using  a simple 
visco-elastic  model  (Reid,  1968) . 

In  the  BEAR  Buoy  program  the  study  of  buoy  motion 
was  used  only  as  a method  for  the  prediction  of  and  sub- 
sequent correction  for  vertical  displacements  of  sensors 
attached  along  the  mooring.  As  such  the  choice  of  problem 
approach  was  tailored  to  the  scientific  objectives  and 
engineering  design.  These  can  be  summarized  as  follows. 
Scientific  objectives ; 

1.  To  study  the  internal  wave  field  using  temperature 
measurements  (thermistors)  as  a means  of  detecting 
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the  internal  wave  signal.  The  thermistors  were 
spaced  along  the  cable  from  near  surface  to  a 
depth  of  approximately  1800  m.  The  emphasis  is 
the  study  of  long  period  internal  waves,  in 
particular,  the  internal  tides. 

Engineering  Design : 

1.  The  individual  thermistor  units  were  potted  and 
spliced  into  the  thermistor  cable.  Each  splice 
was  encased  in  a steel,  egg-shaped  protective 
housing.  The  time  constant  of  the  imbedded 
sensor  was  of  the  order  of  13  min.,  slightly 
greater  than  the  12  min.  sampling  interval. 

2.  The  mooring  sensors  were  mechanically  damped 
to  filter  out  surface  wave  generated  motions 
as  well  as  high  frequency  line  strum  due  to 
vortex  shedding. 

The  above  design  features  were  used  to  minimize 
signal  contamination  due  to  motions  at  frequencies  higher 
than  the  sampling  interval. 

3.  The  mooring  design  was  the  following.  Under 
conditions  of  zero  external  forcing  a non-zero 
tension  was  maintained.  With  increased  external 
forcing  the  dynamic  load  on  the  mooring  was  in- 
creased by  lifting  chain  off  the  bottom.  This 
served  to  reduce  the  magnitude  of  the  horizontal 
excursions  of  the  buoy  and  to  prevent  exotic 
curvatures  in  the  mooring  line. 

Given  the  above  criteria  the  mooring  analysis  was 
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treated  as  a static  problem.  Dynamic  effects  due  to 
higher  frequency  motions  were  assumed  to  be  eliminated 
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by  the  above  engineering  design.  A model  was  developed 
which  utilizes  the  method  of  finite  elements  to  compute 
the  mooring  configuration  in  an  x,z  plane.  The  model 
uses  as  explicit  inputs  prescribed  boundary  conditions 
and  external  forcing  on  the  mooring  system. 

In  application  the  sensor  signal  was  corrected  for 
vertical  displacements  using  a modified  version  of  the 
mooring  configuration  model.  The  modified  version  com- 
putes the  external  forcing  (surface  and  subsurface)  acting 
on  the  system  using  the  measured  quantities  - cc±)le  in- 
clination and  tension  at  the  surface  and  cable  inclina- 
tion at  the  bottom  of  the  thermistor  string  (1825  m)  . 

Given  the  inclination  and  tension  at  the  surface  the 
total  surface  forcing  on  the  buoy  can  be  determined  but 
the  components  due  to  wind,  wave,  and  current  cannot  be 
separated.  In  the  initial  pass,  the  method  uses  the 
surface  input  (inclination  and  tension)  to  compute  the 
surface  force  (total)  and  solves  the  set  of  force- 
balance  equations  assuming  a no-subsurface  drag  ccxidi- 
tion.  If  the  solution  does  not  converge  then  subsurface 
drag  is  added.  Under  a drag  condition  the  horizontal 
tension  is  not  constant  with  depth.  As  such  the  dif- 
ference in  measured  inclination  between  the  surface  eind 
1825  m is  proportional  to  the  sum  of  the  drag  between 
those  depths.  To  satisfy  this  condition  a current  pro- 
file is  computed  which  satisfies  the  drag  requirement. 

A simple  linearly-  decreasing  profile  is  assumed. 


Jd 
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2 . 0 The  Mooring  Configuration  Model  (MCM) 

This  section  treats  the  details  of  the  development 
of  the  mooring  model.  The  mathematical  description  in- 
cludes the  formulation  of  1)  the  surface  forces  of  wind, 
wave,  and  current  acting  on  the  surface  buoy,  2)  the  sub- 
surface drag  on  the  mooring  due  to  current,  and  3)  the 
elasticity  of  the  mooring  components. 

The  model  was  evaluated  using  data  collected  at  the 
time  of  the  buoy  implant  and  during  the  Ground  Truth 
Experiment  (GTE)  conducted  in  April,  1976  (Echternacht , 

23  March  1976) . The  data  included  measurements  of  the 
surface  conditions  (wind,  wave,  and  current),  the  total 
tension  on  the  mooring,  and  XBT  profiles.  The  data 
were  used  to  evaluate  the  predicted  mooring  configura- 
tion by  comparing  the  buoy  temperatures  corrected  for 
vertical  displacements  as  computed  by  the  model  vs  the 
XBT  measurements. 

2 . 1 General  Assumptions  and  Capabilities  of  the  Model 

The  mooring  configuration  model  (MCM)  is  a two 

dimensional  model  used  to  determine  the  configuration  of 
the  mooring  in  an  x-z  plane.  The  model  requires  as  input 
1)  the  buoyancy  force  on  the  buoy,  2)  the  surface  para- 
meter (wind  and  current  speed  and  wave  height)  used  to 
compute  the  surface  forces  acting  on  the  buoy,  and 
3)  a current  profile  through  the  water  column.  As 
discussed  in  Section  1.2  the  mechanical  and  sensor  sub- 
systems were  designed  to  low  pass  signals  with  periods 
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of  the  order  of  minutes  and  longer  thus  eliminating  or 
at  least  minimizing  the  dynamic  motions  of  the  surface 
wave  field  and  surges  in  wind  and  current.  As  such  the 
model  formulation  is  treated  as  a static  problem.  When 
used  to  correct  temperature  measurements  for  signal  con- 
tamination due  to  vertical  displacements  of  the  sensors 
the  mooring  configuration  for  each  data  acquisition  time 
is  solved  separately  (to  be  explained  further  in  Section 
3.0)  . 

An  earth  fixed  Cartesian  coordinate  system  is  used 
with  the  origin  placed  at  the  undisturbed  level  of  the 
free  surface  (buoy  mean  waterline)  and  z is  positive 
downwards.  All  external  forces,  surface  and  subsurface, 
are  assumed  to  act  in  a horizontal  plane  in  either  the 
positive  or  negative  x-direction.  The  mooring  sensors 
(tensiometer  and  inclinometers)  only  measure  the  resultant 
scalar  components  of  the  applied  forces.  No  instrxamenta- 
tion  was  included  to  resolve  individual  force  components; 
i.e,  wind,  wave,  and  current  at  the  surface  and  current 
through  the  water  column.  In  the  model  these  components 
are  added  separately  as  independent  forces.  However,  as 
in  the  case  of  the  horizontal  surface  forces,  it  is  the 
resultcuit  of  these  forces  that  is  used  to  specify  the 
horizontal  component  of  the  surface  boundary  condition. 

The  resultMt  is  expressed  as  a vector  component:  ampli- 
tude and  direction  (positive  or  negative  x) . The  vertical 
component,  to  be  discussed  further  on,  is  the  buoyancy 


force 


The  surface  waves  are  assumed  to  be  small  amplitude, 
deep  water  gravity  waves.  Wave  period  and  subsequent 
displacement  motions  of  the  buoy  due  to  the  wave  field 
are  neglected  (refer  to  Section  1.2).  Hence,  only  the 
amplitude  of  the  wave  field  is  considered. 

All  mooring  lines  are  considered  elastic.  In  the 
formulation  the  nonlinear  elasticity  dependent  on  prior 
loading  history  is  neglected.  A stress-strain  curve  is 
derived  based  on  the  percent  of  load.  Different  load 
curves  are  used  for  the  individual  mooring  components. 

The  load  curves  were  derived  using  data  from  static  load 
tests  conducted  by  the  manufacturers.  Dynamic  effects 
on  the  load  curves  are  neglected.  Internal  damping 
forces  in  the  mooring  lines  are  assumed  negligible  com- 
pared with  viscous  drag  and  stiffness  forces.  Hydro- 
static pressure  forces  are  not  considered.  Sensor 
packages  attached  along  the  mooring  are  treated  as 
added  weight  only.  The  thermistor  units  are  encased 
in  15  protective  cylindrical  egg-shaped  housings.  The 
maximum  housing  diameter  (0.1016  m)  and  length  (0.3483  m) 
compared  to  the  thermistor  cable  dimensions  (0.0252  m 
and  1924  m,  respectively) , is  considered  to  be  a 
negligible  addition  to  the  total  viscous  drag  force. 

Since  the  thermistor  units  are  spaced  at  nearly  regular 
intervals  along  the  cable  the  housing  weights  are  added 
to  the  total  and  evenly  distributed  along  the  cable. 

No  allowance  for  line  creep  of  the  nylon  component  is  in- 
cluded, but  it  is  recognized  that  it  is  possible  over 
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long  periods  of  time.  The  cable  and  nylon  components 
j used  were  torque  balanced  such  that  elongations  due  to 

line  rotation  are  not  considered.  All  wave  forces 
(exciting  and  damping)  acting  on  the  line  are  neglected. 

Viscous  drag  forces  due  to  current  are  computed 
using  the  square  drag  law.  The  forces  are  assumed  to 
act  on  the  line  cross-section  in  the  direction  of  the 

flow.  For  the  buoy  the  wetted  cross-section  is  'used. 

A sectioned  current  profile  is  used.  Each  section  is 

linear  but  most  actual  profile  shapes  can  be  approximated 
by  increasing  the  number  of  sections.  A total  of  n-1 
sections  is  possible;  n is  the  nvunber  of  mooring  elements. 

To  save  computer  time  the  number  of  sections  are  limited 

i 

to  ten  maximum.  j 

2 . .?  The  Surface  Buoy 

The  surface  buoy  has  a total  weight  of  13061  Kg  and 
a positive  buoyancy  of  10715  Kg  per  meter  of  draft.  The 
shape  is  cylindrical  with  dimensions:  3.6  m diameter, 

2.7  m length. 

The  forces  acting  on  the  buoy  are  the  surface 
vertical  force  and  the  horizontal  forces  - viscous  drag 
force  due  to  surface  current,  wind  force,  and  wave 

force.  The  buoy  is  modeled  as  a rigid  surface  piercing  | 

body.  The  coordinate  system  is  shown  schematically  in 

Fig.  2.  The  coordinate  system  used  is  fixed  Cartesian 

with  z positive  downwards.  The  origin  is  placed  at  the 

mean  computed  buoy  water  line.  The  buoy  is  considered 

to  be  upright  at  all  times  (vertical  buoy  axis  parallel 

— — J 
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to  the  z-axis) . The  x-axis  is  by  convection  positive  in 
the  downstream  current  direction. 


Fig.  2 

Coordinate  System  for  the  Buoy 

H is  the  total  length  of  the  buoy,  h©  the  computed  draft, 
Zj^  the  depth  of  penetration,  and  h^  the  buoy  freeboard. 
The  following  expressions  are  used  to  determine  h^  and 
h^,  respectively. 


h = 
o 


Fv  - 


”ab 


B 


h^  = H - h^ 
Si  o 


(1) 


(2) 


where 


Mg, 


is  the  buoy  weight  in  air,  g the  acceleration 
due  to  gravity. 


J 
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is  the  surface  vertical  force  due  to  the 
weight  of  the  mooring,  and 
B the  buoyancy  per  meter  of  draft. 

Eq.  (1)  is  an  expression  for  the  static  equilibrium  posi- 
tion of  the  buoy  as  a function  of  specific  applied  surface 
eind  subsurface  forces.  Knowledge  of  h^  does  not,  however, 
yield  a unique  solution  for  mooring  geometry,  but  rather 
a family  of  solutions  is  possible  depending  on  the  magni- 
tudes of  the  individual  external  forces.  Fig.  3 


Fig.  3 

Condition  of  Static  Equilibrium 


illustrates  schematically  the  static  equilibrium  state. 

F^  is  the  viscous  drag  force  due  to  the  surface  current, 

F the  wind  drag  force  on  the  freeboard  portion  of  the 

a 

buoy,  and  F^  the  wave  force.  is  the  tangential  tension 

at  the  surface  imparted  by  the  mooring  and  9^  is  the 
surface  tension  angle  with  respect  to  the  horizontal. 

The  surface  forces  are  formulated  as  follows. 
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Surface  vertical  force, 

, is  effectively  the  total  weight  of  the  mooring 
acting  on  the  buoy.  For  the  static  case  the  buoy  weight, 
is  balanced  by  the  buoyancy  term,  B (principle  of 
Archimedes).  As  outlined  in  Section  1.2  the  mooring 
system  was  designed  to  provide  a non-zero  tension  under 
conditions  of  zero  surface  and  subsurface  forcing.  The 
details  of  the  design  criteria 'are  provided  in  a latter 
section  (Section  2.6.1).  Under  conditions  of  zero  ex- 
ternal forcing  a minimum  vertical  surface  force  of  4270.0 
Kg  is  maintained  and  the  minimum  draft  is  1.62  m. 

Wind  force,  F^ 

The  wind  force,  F^  was  calculated  using: 

O. 

Fg  = K e*.  A ( Va^j  (3) 

where , 

A is  the  dry  buoy  cross-section, 

V the  wind  speed, 

-3 

0 is  the  density  of  air  (Kg  m ) and 
a 

K the  non-dimensional  drag  coefficient. 

Viscous  drag  force,  F^ 

The  viscous  drag  force  on  the  buoy  and  mooring  line 
follow  a velocity  square  relationship  (Casarella  and  Parsons, 
1970) . In  the  case  of  the  buoy  the  relationship  is  written 


as  follows. 
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F 


d 


<.  A 
2 


(4) 


where , 

is  the  drag  coefficient,  generally  taken 
as  unity, 

mean  sea  water  density, 

A the  wetted  cross-section,  and 

V the  current  speed 


Eq  (4)  defines  the  form  of  the  normal  loading  function. 

In  the  case  of  the  mooring  line  the  normal  loading  relation- 
ship used  is  as  follows. 

C,  J?  DS  sin  e 

(5) 

D is  the  cross-sectional  diameter  of  the  mooring 
line , 

S the  increment  of  line  length,  and 
0 the  inclination  of  the  mooring  line  from 
the  horizontal. 


Tangential  drag  components  are  neglected.  As  shown  by 
Hartman  and  Nash  (1975)  from  test  data  collected  in  the 
Gulf  Stream  tangential  loadings  are  generally  of  the  order 
of  2%  of  the  normal  components. 

Wave  force,  F., 

w 

Wave  forces  are  much  more  complex  and  difficult  to 
model.  The  approach  was  to  represent  the  wave  force  on 
the  buoy  as  the  sum  of  hydrodynamic  drag  and  inertia 
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components  acting  on  a cylinder.  It  is  assumed  that  the 
waves  are  of  small  amplitude,  deep  water  gravity  waves; 
ie,  the  ratio  of  wave  amplitude  to  length  a/L  <<  1,  Thus 
the  wave  force  can  be  written  as  follows  (Dean  and  Hartman, 
1968) . 


D 2 


a cos  <rt  - — 


a sin<j-t 


where , 


Cjj  and  are  the  hydrodynamic  drag  and  inertia 
coefficients,  respectively, 
water  density,  and 
D the  buoy  diameter 


The  above  expression  can  be  simplified  further  by 
eliminating  the  time  dependency.  By  averaging  over  time 
Eg.  (6)  becomes  the  following. 


T I w 


^9CjjD  2 11  2 

(t)dt  = — ^ a ^ \ cos'^^t  dt  - 


4?  gc^n 


if- 


sin  C t dt 


The  first  term  on  the  right  hand  side  is  the  mean  square 
value  and  the  second  integrates  to  zero.  The  wave 
force  is  computed  using  the  following. 


w 4 


2 . 3 The  Mooring  Line 

The  mooring  line  is  comprised  of  three  sections. 


The  upper  section  which  connects  directly  to  the  buoy 
is  the  1-inch  diameter  steel  thermistor  cable.  The 
center  section  is  a 2-  inch  nylon  mooring  line  and  the 
bottom  section  1-1/8  inch  stud  link  chain.  It  should 
be  noted  that  under  conditions  of  zero  external  forcing 
approximately  42  m of  chain  are  lifted  off  the  bottom. 
During  conditions  of  non-zero  external  forcing  the  static 
load  increases  due  to  the  lifting  of  additional  chain  off 
the  bottom.  The  model  foirmulation  includes  the  changes  in 
static  load  and  the  subsequent  changes  in  elongation  of  the 
nylon  and  cable  members  due  to  varying  conditions  of  ex- 
ternal forcing. 

The  mathematical  model  of  the  statics  of  the  mooring 
line  is  formulated  using  the  method  of  finite  elements. 

In  that  approach,  the  mooring  line  differential  equations 
are  integrated  incrementally  down  the  mooring  line  to  ob- 
tain the  static  equilibrium  position  for  each  line  element. 
The  spatial  coordinate  of  each  element  is  computed  with 
respect  to  the  origin  located  at  the  buoy.  The  viscous 
drag  force  on  each  line  elemement  due  to  the  velocity  field 
is  included. 

In  the  finite  element  formulation  the  mooring  line  is 
reduced  to  a number  of  s-raight  line  segments  joined  at  end 
points  or  nodes.  Each  segment  is  considered  an  extensible 
cable  without  mass;  i.e.. 


a straight  spring.  All  forces 
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J 

i 

■i 


and  mass  along  the  cable  are  assumed  to  be  concentrated  at 
the  nodes.  Fig.  4 shows  schematically  the  representation 
of  the  mooring  line. 


Finite  Element  Formulation  of  the 
Mooring  Line 

To  illustrate  the  formulation,  consider  the  balance 
of  forces  acting  on  a differential  element  of  the  mooring. 
As  shown  in  Fig.  5 the  line  forces  are:  (1)  the  constant 
force  due  to  gravitation  and  (2)  the  variable  tension 
forces  transmitted  from  adjoining  elements,  and 
The  constant  force  due  to  gravitational  attraction  is  the 
wot  weight  per  unit  stretched  length  of  element  i.  The 
fluid  force  is  the  viscous  drag  force,  F^^,  duo  to  the 
fluid  velocity  as  given  by  Eg.  (5)  (Section  2.2).  For 
static  equilibrium: 


Ti  t . (mg),  t t^^^ 


j 


0 


(9) 
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Fig . 5 

Forces  Acting  on  Line  Mass  Element  (i) 

Eq.  (9)  is  solved  sequentially  for  each  mass  element  i; 
i = 1 at  the  surface,  n at  the  bottom.  Near  the  surface 
the  dominant  component  of  the  tension  is  the  vertical  com- 
ponent due  to  the  v?eight  of  the  mooring  below  element  i. 

With  increasing  depth  the  vertical  component  decreases  and 
the  horizontal  component  due  to  drag  increases.  Near  the 
bottom  the  horizontal  component  dominates  resulting  in  near 
horizontal  line  angles. 

The  horizontal  component  due  to  drag,  as  given  by  Eq. 

(5) , is  computed  using  V^.  Here  is  the  average  current 

magnitude  acting  on  the  differential  length  of  line  element, 

S^.  is  taken  from  the  current  profile. 

The  tensile  forces  acting  on  the  n mass  elements  and 

the  inclination  of  the  line  elements  change  with  depth  to 

keep  all  the  above  forces  in  equilibrium  (refer  to  Fig.  5). 

Consider  the  following  example.  Given  the  equilibrium  posi- 
A 

tion,  Y^,  of  the  ith  element,  then  the  inclination  of  the 

tension  vector,  acting  on  the  ith  element  and  the 

A 

equilibrium  position,  of  mass  element  i+1  can  be 
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found  as  follows.  Y is  the  geometric  position  vector. 

Since  the  vertical  component  of  the  tension  is  simply 
the  weight  of  that  portion  of  the  mooring  acting  on  mass 
element  i,  then 


F ...  = sin  9 

v,i+l  s s 


j. 

-/ 

j=l 


(mg)j  dz 


and 


cos  9^  + 

h,  1+1  s s 


1 

I 

j=l 


Fdj  dz 


where , 


(10) 


F • and  F,  . are  the  vertical  and  horizontal 
V , i+i  h , 1+1 

components  of  the  tension,  and  the  tension  at  the 

surface  (buoy) . 

Thus  the  tension,  inclination  angle,  and  equilibrium 
position  of  element  i+1  are  as  follows. 


i+1  ~ {^v,i 


2 

i+1 


+ F 


h,ii-l/ 


(11) 


and 


®i+l  ~ ^^v,i+l/\,i+l^ 


X,  i+1 

Y 

X,1 

^ ^i+1 

cos 

®i+l 

z,i+l 

Y 

Z,1 

®i+l 

sin 

®i+l 

(12) 


(13) 


Here  x,  z refer  to  the  x,z  component  directions  amd 
is  the  unstretched  length  of  line  element  i+1.  It 
should  be  pointed  out  at  this  time  that  line  stretch  is 
included  in  the  formulation  used  to  determined  the  geometric 
position.  The  details  of  the  formulation  will  be  given  in 
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the  next  section. 

To  briefly  summarize,  the  mooring  configuration  is 
determined  as  follows.  The  force  balance  equation  is  com- 
puted sequentially  for  line  elements  i=l,  n to  determine 
the  equilibrium  position  of  each  mass  element  thus  describing 
the  geometrical  shape  of  the  mooring  in  an  x,z  plane.  The 
advantage  of  this  method  is  that  the  equations  can  be  derived 
directly  from  Newton's  Laws  of  motion.  However,  it  should 
be  noted  that  this  method  is  limited  to  static  problems  with 
prescribed  boundary  conditions  and  external  forces. 


2 . 4 Mooring  Elasticity 

To  accurately  predict  the  mooring  configuration  the  line 
stretch  must  be  included.  The  elasticity  of  the  mooring  for 
the  cable  and  nylon  sections  under  varying  load  was  included 
explicitly  in  the  model  formulation.  It  was  assumed  that 
the  materials  followed  Hooke's  Law;  ie.,  the  extensibility 
of  the  mooring  line  is  proportional  to  the  pulling  on  the 
line  due  to  tension.  Extensibility  by  squeezing  due  to 
hydrostatic  pressure  was  not  considered. 

Having  determined  the  tension  acting  on  the  i+1  element 
of  line  (Eq.  11)  the  elastic  stretch  is  calculated  as  follows. 
From  Hooke's  Law  the  incremental  stretch,  ,6  is  defined  as: 

'W 


Xn  is  the  modulus  of  elasticity  with  dimensions  of  force  per 
unit  cross-sectional  area.  A;  ie,  /Ng  = A 'A  . The  modulus  X 
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was  determined  from  static  bench  tests  conducted  by  the 
manufacturers.  The  manufacturers'  test  data  are  given 


in  Table  2. 

For  the  nylon 

, eight  strand 

pli-moor 

Double 

Armored  Cable 

% 

Tension 

Modulus , 

Stretch 

% 

Load 

(Kg) 

Ac 

(m) 

Stretch 

3.5 

1104.1 

4.5x10^ 

6.0 

.31 

6.7 

2268.0 

8.1x10^ 

6.8 

.36 

20.0 

6804.0 

1.2x10® 

13.7 

.71 

33.6 

11339.9 

1.6x10® 

17.2 

.89 

43.0 

14601.0 

1.7x10® 

21.0 

1.09 

Eight 

Strand  Pli-moor 

Nylon 

% 

Load 

Tension 

(Kg) 

Modulus , 

AN 

Stretch 

(m) 

% 

Stretch 

2.0 

834.6 

7.59x10^ 

99.4 

3.66 

5.0 

2086.5 

1.02x10^ 

172.8 

6.37 

10.0 

4173.1 

1.23x10^^ 

294.1 

10.84 

15.0 

6259.6 

1.37x10^ 

397.1 

14.64 

40.0 

16692.4 

2.25x10^ 

640.9 

23.63 

Table  2 

Memufacturer  Bench  Test  Data  - Load  vs 
Percent  Stretch:  Cable  and  Nylon 
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construction,  X varies  from  7.6  x 10^  at  2%  loading  to  greater 

4 

than  2.2  X 10  at  leading  greater  than  40%.  For  the  double 
armored  conductor  cable  X varies  from  4.5  x 10^  at  3.5% 
loading  to  greater  than  1.6  x 10^  at  loads  greater  than  43%. 

Values  for  the  modulus  X for  both  the  cable  ( Xc)  and 
nylon  (Xj^)  members  were  computed  from  the  test  data  in  one 
percent  intervals  over  the  rainges  given  in  Table  2.  The  values 
for  .Xjj  cind  Xc  are  shown  in  Fig.  6.  In  the  model  X was  computed 
as  follows,  (1)  Percent  load  was  computed  for  each  line  element 
i using: 

% Load  = (T/LC)  • 100  (15) 

where 

LC,  the  load  constant,  is  the  rated  breaking  strength  - 
3.385  X 10^  Kg  for  the  cable  cmd  4.173  x 10^  Kg 
for  the  nylon . 

(2)  the  appropriate  X corresponding  to  Eq.  (15)  was  found  using 
a table  look  - up. 

Given  X the  incremental  stretch,  was  computed 

using  Eq.  (14).  Once/iS^^j^,  was  determined  the  element  co- 
ordinates (given  by  Eq.  13)  were  corrected  to  include  the 

stretch. 

/ 

Y 

x,i+l 
^z,i+l 


J 


(16) 


The  prime  denotes  the  stretch  included  coordinate. 


TENSION  (Kg  X 10*) 


Load  and  Modulus  Curves 
Nylon  and  Cable 


MODULUS.  X 


wmfmmm 
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2 . 5 Method  of  Solution 

The  previous  sections  outlined  the  necessary  ingredients 
for  the  mathematical  model.  This  section  details  their  useage 
in  obtaining  a solution  for  given  external  inputs. 

The  model  requires  as  input  the  buoyant  force  on  the 
buoy;  surface  environmental  parameters  of  wind  speed,  wave 
height,  and  surface  current  velocity;  cind  a current  profile 
through  the  water  column.  The  environmental  surface  para- 
meters are  used  to  compute  the  surface  horizcsital  force  com- 
ponent acting  on  the  buoy.  For  the  static  solution,  the 
analysis  is  started  at  the  buoy.  The  surface  force  com- 
ponents are  used  to  compute  the  line  tension  and  angle  of 
inclination  at  buoy.  These  are  the  imposed  surface  boundary 
conditions  and  act  at  the  junction  between  the  buoy  cind  the 
first  element  of  the  n element  elastic  mooring  line.  The 
pertinent  static  equations  for  the  buoy  (Section  2.2)  with 
the  attached  mooring  line  (Sections  2.3  and  2.4)  are  as 
follows. 

F - T = 0 (17) 

X X 

''v  - ”g  - ^2  ' “ 


where 


F 


X 


M 


g 


is  the  sum  of  the  surface  forces, 
F,  + F,  + F,^,  and 

Q ci  W 

is  the  weight  of  the  mooring. 
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By  specifying  we  have  three  unknowns:  F^,  and  T^. 

Initially,  we  assume  a value  for  the  surface  vertical  force, 
/ 

F^,  Where, 


Thus  the  unknowns  can  be  solved  using  the  following, 

0g  = tan'^ 


2 + F'2)^i 


T = (F  + F 

S X V 


T = sin 
z s s 


The  subscript  s denotes  surface  parameters. 

Having  solved  the  buoy  force  the  mooring  line  is  analyzed 
next.  For  the  mooring  line  formulation,  a numerical  integra- 
tion is  carried  down  the  line  in  a set  of  n finite  segments. 

In  this  representation  the  line  equations  can  be  written  as 
follow;3. 

'■x.i  - ’■x.l  - ''di  ' “ 

"v,!  - - ’■z.i  ' “ 

F^^  is  the  current  drag  on  element  i computed  from  the 
current  profile. 

Solving  the  force  balance  equations,  (23)  and  (24),  for 
element  i the  conditions  for  the  i+1  element  are  specified 
as  follows. 

’■i+l  - '''v.i'  * '’x.i^'’’  <“> 
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The  coordinate  position  of  each  element  Is  found  using 
Eqs.  (13)  and  (16).  The  iteration  is  carried  out  to  the 
end  of  the  mooring;  1 ■>  n.  At  that  point  the  bottom 
boundeiry  conditions  should  be  met,  such  that: 


F 


v,n 


0 


and 


(27) 


n 


i=l 

Z'  is  the  measured  depth  = 4793  m,  and 

£ = 0.001. 

Thus  if  the  computed  depth  does  not  agree  to  within  0.1% 
the  process  is  repeated  with  a new  F^'  specified. 


2.6  Model  Case  Studies 

The  model  (MCM)  was  used  to  study  the  following.  (1) 

Implant  simulation  - studies  were  conducted  for  conditions 

of  zero-external  forcing  over  a ramge  of  water  depths  and 

\ 

varying  mooring  line  lengths.  These  studies  were  conducted 
prior  to  the  buoy  implemt  to  determine  the  optimum  mooring 
length  for  the  cuiticipated  range  of  bottom  depths.  (2)  En- 
vironmental simulation  - these  studies  were  run  to  determine 
the  mooring  configuration  under  Veirying  conditions  of  ex- 
ternal forcing  by  the  environment. 

2.6.1  System  Implant  Simulation 

Prior  to  implant  a bottom  survey  was  conducted  to  deter- 
mine the  depth  of  bottom  and  chemges  in  bottom  topography 
within  the  proposed  study  area.  The  Implant  site  was 
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located  on  the  abyssal  plain  approximately  15  Km  seaward  of 
the  continental  margin  at  a depth  of  approximately  4800  m. 
The  bottom  is  relatively  flat  with  minor  topographic  varia- 
tions - of  the  order  of  a few  meters.  The  mooring  design 
was  intended  as  a semitaut  system  with  a non-zero  minimum 
tension  (refer  to  Section  1.2).  To  achieve  this  the  design 
specified  that  a short  section  (of  the  order  of  tens  of 
meters)  of  ground  tackle  (chain)  would  be  maintained  off 
the  bottom  under  conditions  of  zero  - external  forcing. 

As  such  a knowledge  of  the  mooring  configuration  under 
these  conditions  was  necessary  in  order  to  specify  the 
lengths  of  the  mooring  members  (cable  and  nylon) . 

Using  the  MCM  a number  of  case  studies  were  run. 

The  cases  included  (1)  varying  the  bottom  depth  with 
mooring  lengths  held  constant  and  (2)  varying  mooring 
lengths  with  depth  held  constant.  The  results  of  all 
case  studies  will  not  be  presented  only  the  case  repre- 
senting the  final  mooring  specification  and  measured  depth 
at  the  implant  site. 

The  input  specifications  for  the  final  case  were  the 
following:  (1)  bottom  depth  - 4800  m,  (2)  cable  length  - 
1924.4  m,  and  (3)  nylon  length  - 2713.6  m.  Fig.  7 illus- 
trates the  results.  In  the  figure,  = 0 represents  the 
depth  - 4800  m.  For  this  implant  depth  approximately  42  m 
of  chain  are  maintained  off  the  bottom  with  a minimum  non- 
zero tension  of  approximately  4270  Kg.  It  should  be  noted 
that  stretch  in  the  nylon  and  cable  members  was  included  as 
shown  in  the  figure.  The  final  mooring  lengths  as  specified 
by  this  case  study  allowed  a safety  margin  of  approximately 
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145  m in  the  depth  at  the  implant  site.  The  results  of 
the  case  study  and  thus  the  MCM  calculations  were  verified 
at  a later  date  by  measuring  the  weight  of  the  mooring. 

During  the  Ground  Truth  Experiment  conducted  in  April, 

1976  (Echternacht , 23  March  1976)  the  mooring  weight  was 
measured.  The  sea  and  weather  states  were  calm.  The  method 
entailed  hauling  on  the  mooring  until  the  buoy  tensiometer 
was  zero  and  maximum  on  the  hauling  tensiometer.  The  mooring 
weight  - measured  tension  - agreed  with  the  total  calculated 
model  line  weight  - thus  line  length  - within  the  error 
limits  of  the  instrumentation  used. 

2.6.2  Environmental  Simulation 

This  section  presents  the  results  of  specific  case 
studies  using  the  MCM.  The  case  studies  to  be  presented 
represent  a wide  range  of  external  forcing  by  the  environ- 
ment. The  cases  presented  are  separated  into  three  categories: 
surface  forcing  in  the  presence  of  a current  - linearly  de- 
creasing with  depth  and  with  shear  - and  surface  forcing  in 
the  absence  of  a current.  The  above  categories  were  used, 
in  part,  to  assess  the  effect  of  subsurface  drag  on  the 
mooring.  The  application  of  those  results  are  contained  in 
a latter  section  which  details  the  formulation  of  the  model 
used  to  correct  sensors  (thermistors)  for  vertical  displace- 
ments . 

Table  3 lists  the  three  case  study  categories  and  the 
magnitudes  of  the  environmental  inputs.  Within  each  category 
the  inputs  represent  conditions  of  strong,  moderate,  and  mild 
forcing  covering  the  range  of  expected  sea  and  weather  states 


Case 

Study 

Wind 

(nis“^) 

Wave 

(m) 

Current 
sfc  (ms"^) 

Current 

Profile 

1 

10.30 

3.05 

1.03 

1.03  - 0. 

2 

6.90 

2.00 

.66 

.66  - 0. 

3 

3.50 

1.00 

.33 

.33  - 0. 

4 

10.30 

3.05 

0. 

0. 

5 

6.90 

2.00 

0. 

0. 

6 

3.50 

1.00 

0. 

0. 

7 

5.00 

2.00 

1.03 

1.03-(-.5)  (2500) 
(-.5)-0.  (4793) 

8 

5.00 

2.00 

.50 

( .5)-(-.5) (2500 
(-.5)-0.  (4793) 

9 

5.00 

2.00 

Table 

0. 

3 

0.-(-.5)  (2500) 

(-.5)-0.  (4793) 
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with  the  exception  of  rare  events  such  as  hurricanes. 

Cases  1-3  and  7-9  differ  in  the  typo  of  current  profile 
used.  The  former  uses  a linearly  decreasing  profile  with 
surface  maximum;  the  latter  includes  a current  shear 
(variable)  under  conditions  of  moderate  surface  forcing. 

In  order  to  evaluate  the  effect  of  external  forcing 
for  the  various  cases  the  vortical  displacement  of  reference 
sensors  within  specified  zones  was  used.  The  depth  zones, 
as  given  in  Table  4,  correspond  to  regions  of  nearly  constant 
vertical  temperature  gradient,  The  zones  were  taken 

from  an  average  temperature  profile  based  on  XBT  data  collected 
by  the  U/V  I.AMB  in  the  study  area  during  the  early  phase  of 
the  experiment.  The  vertical  displacement  is  computed  witli 
respect  to  a reference  or  nominal  depth.  Here  the  nominal 
depth  is  defined  as  tl^e  reference  sensor  depth  for  a specific 
boundary  condition  of  surface  tension  excluding  horizontal 
displacements  of  the  mooring.  Thus  the  nominal  depth  in- 
cludes changes  in  mooring  length  due  to  stretch  resulting 
from  changes  in  total  mooring  weight.  As  an  example,  typical 

values  of  surface  tension  for  conditions  of  strong  and 

4 3 

moderate  forcing  are  10  and  7x10  Kg,  respectively.  The 
corresponding  nominal  depths  for  the  zone  6 reference 
sensor  arc  1854.5  and  1843.2  m.  For  the  zero-external 
forcing  state  the  nominal  depth  for  the  zone  6 reference 
is  1832.1  m.  Thus  under  strong  forcing  the  reference  sensor 
is  22.1  m deeper  duo  to  cable  stretch.  The  vertical  excursions, 
^z,  from  nominal  result  from  changes  in  mooring  configuration 
when  horizontal  forcing  is  included. 


M 


i 
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The  results  of  the  simulation  studies  corresponding  to 
the  external  inputs  given  in  Table  3 are  listed  by  category 
in  Table  5A  through  C.  For  the  current  shear  cases  (7-9)  the 
depth  at  the  inflection  point  is  in  parenthesis. 


CASE 


Depth 


Surface  Forcing  - Subsurface  Drag 


CASE 


Surface,  Subsurface  Forcing  - with  current  shear 
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3.0  The  Vertical  Displacement  Corrector 


This  section  treats  the  methodology  used  to  correct 
the  signal  of  sensors  attached  along  the  mooring  line  for 
contamination  resulting  from  changes  in  mooring  configura- 
tion. The  first  two  parts  include  1)  the  detailed  develop- 
ment of  the  vertical  displacement  corrector  model  (VDCM) 
and  2)  the  methodology  used  to  compute  the  vertical  dis- 
placements of  the  individual  sensors  and,  subsequent!-' , 
correction  of  the  sensor  signal.  The  final  part  of  this 
section  examines  the  model  application.  Included  is 
1)  a comparison  of  the  configuration  computed  - VDCM  vs 
MCM  - for  varying  boundary  conditions  and  2)  a comparison 
of  the  corrected  signal  vs  XBT  measurements. 


3 . 1 General  Assumptions  and  Capabilities 

The  VDCM  is  a modified  version  of  the  mooring  con- 
figuration model.  The  assumptions  used  to  mathematically 
describe  the  mooring  system  are  the  same  as  those  used  for 
the  MCM.  The  method  of  solution  utilizing  the  finite 
element  approach  is  identical  as  well  with  the  exception  of 
the  iteration  scheme  used  to  close  the  solution.  The  primary 
difference  between  the  two  models  is  in  the  approach  used  to 
specify  the  boundary  conditions.  To  recap,  the  MCM  uses 
explicit  inputs  of  1)  the  magnitude  of  the  surface  wind, 
wave,  and  current  to  compute  the  surface  forces  acting  on 
the  buoy  and  2)  the  current  profile  used  to  compute  the 
subsurface  drag.  As  such  all  external  forces  are  given. 

For  the  VDCM,  however,  the  explicit  inputs  listed 
above  are  not  available  and  mv^st  be  derived  using  measured 
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quantities  - the  cable  inclination  and  tension  at  the 
surface  and  the  cable  inclination  at  the  junction  between 
the  cable  and  nylon  members  (approximately  at  a nominal 
depth  of  1900  m) . Given  the  inclination  and  tension  at 
the  surface  the  total  surface  forcing  on  the  buoy  can  be 
determined  but  the  components  due  to  wind,  wave,  and 
current  cannot  be  separated.  As  such,  in  the  initial 
pass  the  method  uses  the  surface  input  (inclination  and 
tension)  to  compute  the  surface  force  (total)  and  solves 
the  set  of  force-balance  equations  assuming  a no-subsurface 
drag  condition.  It  should  be  noted  that  the  solution  is 
unique  for  a given  value  of  the  surface  boundary  condi- 
tions and  the  solution  will  converge  if  and  only  if  the 
no-drag  condition  is  satisfied.  Mathematically  the 
formulation  is  similar  to  the  catenary  problem.  Under 
conditions  of  zero-subsurface  forcing  (here  defined  as 
the  zeroth  state)  the  horizontal  tension  is  constant  with 
depth.  It  should  be  pointed  out,  however,  that  the  in- 
clinations at  the  surface  and  1900  m need  not  be  equal 
to  satisfy  the  zeroth  condition.  If  the  solution  does 
not  converge  then  subsurface  drag  must  be  added.  Under 
a drag  condition  the  horizontal  tension  is  not  constant 
with  depth.  Thus  the  difference  in  measured  inclination 
between  the  surface  and  1900  m is  proportional  to  the 

sum  of  the  drag  between  those  depths  and  the  drag,  in 

2 

turn,  is  proportional  to  v (horizontal  current  speed) . 

To  satisfy  this  condition  a current  profile  is  computed 
which  satisfies  the  drag  requirement.  A simple  linearly- 
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decreasing  profile  is  assumed. 

Formulation  of  Subsurface  Drag  and  Current  Profile 

For  cases  in  which  Lhe  zeroth  condition  is  not  satisfied 
ie.,  the  solution  does  not  converge  - this  implies  that  an 
unknown  percentage  of  the  horizontal  component  of  tension 
is  due  to  drag.  The  actual  current  profile  is  unknown  and 
therefore  the  distribution  of  drag  with  depth.  The  static 
formulation  is  similar  to  the  catenary  problem  with  the 
exception  that  under  drag  conditions  the  horizontal  com- 
ponent of  tension  for  an  arbitrary  depth  is 


i 


j = l 


rather  than  a constant  Fj^  with  depth.  In  Eq.  (28)  the 
subscript  s denotes  the  surface.  To  approximate  the 
integrated  drag  portion  of  Eq.  (28)  the  following  method 
is  used. 

Using  the  measured  line  inclinations  at  the  surface, 

0 , and  bottom  of  the  thermistor  cable,  0 , the  sum  of 
s o 

the  drag  between  the  surface  and  cable  bottom,  by  re- 
writing Eq.  (28)  becomes 


j = l 

For  the  static  case  the  vertical  component  of  tension  at 
Z^  is  simply  the  weight  of  the  mooring  below  Z^.  Thus, 
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I 

I 


F = F 
v,c  V 


(mg)^  dz 


j = l 


and  the  tangential  tension  is: 


T = F (sin  9 ) 
c V,  c c 


-1 


Therefore,  it  follows  that 


F.  = T cos  9„ 
h,c  c c 


(30) 


(31) 


(32) 


As  an  aside  it  should  be  noted  at  this  point  that  the 
solution  of  Eg.  (29)  provides  not  only  the  magnitude  of 
the  sum  of  the  drag  to  Z^,  but  also  the  sign  of  the 
integrated  current.  Thus,  the  following  criteria: 


V 

> 0, 

v > 0 

J 

= 0, 

zeroth  case 

j=l 

A 

o 

V < 0 

As  an  example,  for  v < 0 the  current  increases  the 
horizontal  component  of  tension  in  the  negative  sense 
with  respect  to  the  surface  value.  The  effect  is  to 
reduce  the  magnitude  of  the  horizontal  excursion  of  the 
mooring  and,  in  the  case  of  a strong  current,  can  change 
the  direction  of  the  excursion. 

At  this  point  no  further  assumptions  have  been  applied. 


Since  current  profiles  were  not  measured  the  variability  of 
the  drag  with  depth  is  unknown.  Therefore,  to  simplify 
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the  problem  the  current  profile  is  assumed  to  be  linearly 
decreasing  with  depth  from  a surface  maximum  to  zerj  at 
the  ocean  bottom. 

To  reiterate,  the  drag  acting  on  line  segment  i is 
given  by  Eq.  (5)  - refer  to  Section  2.2. 


c, D D S.  V.  a 

„ _ d^  1 1 sin  9. 

^di 2 ^ 

Thus,  Eq.  (29)  can  be  expressed  as 


^dj 


^d5-  ^ 


C S.v- 

J 3 J 

j = l 


sin  9j  dz 


c^  and  line  cross-section  D are  constant  and  for  simplicity 
changes  of  density  with  depth  are  ignored.  Frcxn  Section 
2.5  it  was  shown  that  for  the  cable  at  43%  load  the  percent 
stretch  is  1.09%  (Table  2).  As  such,  changes  in 
cable  length  have  a negligible  effect  on  the  numerical 
value  of  Eq.  (34) . To  simplify  further,  since  a linear 
current  profile  is  assumed  the  change  in  inclination,  9^, 
is  assumed  to  change  linearly  from  the  surface  to  Z^. 

Thus  Eq.  (34)  becomes 


2i\ 

dz-5- 


where,  we  let 


K denote  the  numerical  difference  in  horizontal 
o 

tensions  between  the  surface  and  , and 
_ c , 0 D S ^ 

^o  “ -S-2 sin  9 (36) 

9 is  the  numerical  average  taken  frcm  the  linearly 
changing  9^  profile  and  S the  line  length. 
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To  construct  the  linear  current  profile  the  average 
value,  V,  of  the  current  between  the  surface  and  is 
at  depth  Z which  lies  midway  between  z = 0 and  Z^.  To 
find  V we  assume  that 


f(  ■ 

) 3 


— I ■ -1 

7 = !^ 


where 

n^  is  the  number  of  segments  between  Z = 0,  Z^. 

The  above  expression  is  not,  of  course,  exact  but  is  a 
close  approximation  for  n^,  large.  In  the  VDCM  n^  is 
generally  0?374. 

Having  determined  V,  the  slope,  m^,  of  the  current 
profile  between  Z and  the  ocean  bottom,  Z^  is 


c Zb  - Z 

and  the  magnitude  of  the  surface  current  is 


V = V + m Z (39) 

s c 

The  sign  of  is  determined  using  the  criteria  given  by 
Eq.  (33) . The  VDCM  is  then  run  with  drag  computed  using 
the  estimated  current  profile. 


3.2  Method  of  Solution 

The  method  of  solution  for  the  VDCM  is  the  same  as 
utilized  by  the  MCM  with  the  exceptions  as  developed  in 
the  previous  section.  With  the  estimated  current  profile 
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I; 


determined  the  models  are  the  same  except  for  the  itera- 
tion method  to  be  given  later  in  this  section.  Eqs  (17) 
and  (18)  are  used  to  determine  the  forces  acting  on  the 
buoy  and  Eqs  (23)  and  (24)  for  the  mooring  line. 

The  iteration  method  is  the  following.  At  the  end 
of  each  iteration  the  bottom  boundary  conditions  are 
applied  in  the  seime  manner  as  was  done  with  the  MCM; 
ie,  Eg.  (27)  - Section  2.5.  For  the  initial  iteration 
the  zeroth  state  is  assumed.  If  the  zeroth  criteria 
are  not  met  the  second  iteration  uses  the  estimated 
current  profile  as  given  in  Section  3.1.  Again  if  the 
computed  depth  does  not  agree  to  within  0.1%  the  process 
is  repeated.  For  the  MCM  the  iteration  variable  is 
but  for  the  VDCM  it  is  Vs  and  the  corresponding  ccanputed 
current  profile.  As  derived  in  Section  3.1  is  a 
function  of  V,  Eg.  (37)  , which  is  an  approximation.  The 
force  balance  formulation  is  sensitive  to  small  changes 
in  horizontal  forcing  - eg,  Eq.  (23)  - such  that  Eq.  (37) 
is  used  only  for  the  initial  estimate  on  the  second  itera- 
tion. Initial  values,  V^,  used  in  subsequent  iterations 
are  found  by  fine  - tuning  the  initial  via  extrapolation 
using  the  following  criteria. 

For,  (Zp  - Z')„>  I^Z'I  and 


(Z„  - z') 

p m 

and 

< 

"p  - 

2’)m-l/^''s 

< 0 

■ Z')„ 
p m 

<“ 

and 

> 

- 

1 

< 

(0 

> 0 

(40) 
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where, 

Z is  the  predicted  bottom  depth, 

P 

Z'  is  the  measured  depth,  and  subscript  m 
the  iteration  number. 

The  new  initial  guess,  V^',  is  specified  by 
'"s’  = ""s  + ^^s 

In  application  the  VDCM  is  run  for  each  acquisition 
period.  The  model  outputs  the  mooring  configuration  and 
the  computed  depth  at  each  sensor  location. 

Computing  Vertical  Displacements  and  Signal  Correction 
To  compute  the  vertical  displacement  at  each  sensor 
location  the  computed  location  depth  as  returned  by  the 
VDCM  was  substr acted  from  the  computed  nominal  depth. 

As  defined  in  Section  2.6.2  the  nominal  depth  is  the 
reference  at  which  the  sensor  would  reside  in  the 
absence  of  changes  in  mooring  configuration.  As  such 
the  sensor  signal  (temperature)  contained  fluctuations 
due  both  to  the  environment  and  to  changes  in  depth. 

To  correct  the  signal  for  bias  due  to  depth  changes  the 
following  method  was  used. 

(1)  The  change  in  temperature  due  to  depth  changes, 

A 

^I’i  = ^T  • ^^i  (42) 

is  the  average  temperature  gradient  for  the 
experimental  area  based  on  historical  data 


DEPTH  (m) 
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(FNWC  viml  XBT  from  the  R/V  LAMB)  . The 
uaod  is  shown  in  Fiq.  8. 

is  the  computed  depth  chaiuie  from  nominal 

.It  sensor  location  i. 

(2)  In  situ  ti'mperature , T' 

T'  = T - ar.  (43) 

m " i 

T is  the  measured  temperature  not  corrected  for 
m 

d i s p 1 a ceme  n t bias. 

3 . 3 Mod^  Ver  i f icat  iott  ^ Corrector  vs  MC.M  - vs  XBT  Pat  a 

Ideally,  verification  of  the  mcxiel  would  include  a 
comparison  of  imxiel  predictions  vs  environmental 
measurements  over  a wide  ranqe  of  environmental  conditions. 
However,  for  the  remote  BKAU  Buoy  system  the  above 
approach  was  not  firactical.  As  an  alternative,  several 
case  studies  were  conducted  to  provide  environmental 
data  which  could  be  used  to  evaluate  the  followinq. 

1)  To  verify  the  capability  of  the  model  to  accurately 
define  the  nominal  zero-external  forcing  state. 

2)  To  verify  tlie  accuracy  cf  the  nuxlel  in  predicting 
mooring  configuration  during  perioiis  of  non-zero 
external  forcing. 

The  model  verification  included  comparisons  of  MCM  vs 
tensiometer  me.isurement  s , MCM  vs  VOCM  results,  and  VDCM 
vs  XBT  measurements. 

Nominal  Cond i t ion 

As  discusseil  in  Section  2.0.1  the  mooring  weight 
was  measured  during  the  C, round  Trutli  Fxperiment.  For  the 
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condition  of  zero-external  forcing  the  measured  surface 
tension  is  equal  to  the  weight  of  the  mooring  which 
includes  the  cable  and  nylon  members  and  a short  length 
of  the  ground  tackle.  As  such,  the  measured  weight 
yields  the  total  mooring  length  with  stretch  included 
since  the  weight  per  unit  length  of  each  member  is 
known.  The  MCM  was  run  using  the  measured  tension  as 
input.  The  model  prediction  of  total  mooring  length 
(stretch  included)  agreed  with  the  field  measurements 
within  the  error  limits  of  the  instrumentation  used. 
Conclusion:  the  model  accurately  predicts  mooring  length 
with  elasticity  included  for  the  nominal  state. 

Externa 1-Forcing  Condition 

To  test  the  accuracy  of  the  model  in  predicting  the 
mooring  configuration  during  the  non-zero  external  forcing 
state,  two  study  comparisons  were  made. 

1)  MCM  vs  VDCM  comparison 

2)  VDCM  vs  XBT  measurements 

MCM  vs  VDCM  Comparison 

The  MCM  and  VDCM  were  compared  for  several  case  studies 
of  differing  conditions  of  surface  and  subsurface  forcing. 
Table  3 (Section  2.6.2)  summarizes  the  inputs  for  each  case. 
For  the  MCM/VDCM  comparison  cases  7-9  were  not  run.  The 


Depth]  AZ  lAT  ATE  AT-  ^ TE  AZ  j AT  ATE 

Zone  (m)  | (Oc)  {°cf  (m)  (Oc)  (Oc7  (Oc)  ,o 


Surface  and  Subsurface  Forcing 
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thermistor  displacement  (referenced  to  the  zero  forcing  - 
nominal  state)  for  each  sensor  location  was  computed  by 
both  the  MCM  and  VDCM  and  compared  for  each  case  study. 
Tables  6A  and  B summarize,  respectively,  the  results  for 
surface  plus  subsurface  forcing  and  for  surface  forcing  - 
zero  subsurface  drag  cases.  For  the  sake  of  brevity, 
only  one  sensor  location  was  chosen  for  each  of  six 
depth  zones  defined  in  Table  4 (Section  2.6.2).  In 
Tables  6A  and  B,^Z  represents  the  vertical  displacement 
(m)  of  the  thermistor  from  the  nominal  (zero  force)  depth 
position.  The  AT  represents  the  temperature  differential 
due  to  vertical  displacement  - referenced  to  the  MCM 
prediction  and  ATE^  is  the  difference  in  differential 
temperatures  between  the  MCM  and  VDCM.  As  seen  in  Table 
6B  the  error  is  essentially  zero  since  there  is  no 
ambiguity  between  the  two  models  for  the  zero  subsurface 
drag  case. 

In  the  comparison  given,  the  depth  zones  (Table  4) 
represent  regions  of  average  gradient.  In  actual 
practice  the  gradients  are  obtained  from  the  average 
temperature  profile  over  a much  smaller  vertical  distance  - 
generally  of  the  order  of  the  nominal  depth  (thermistor 
N)  +100  m,  maximum.  In  comparing  these  values  with  ATE^, 
it  should  be  noted  that  the  error  between  the  MCM  and 
VDCM  is  less  that  the  resolution  (°c/bit)  of  the  thermistors 


used. 


VDCM  VS  XBT  Measurements 

The  model  results  using  acquired  data  were  compared 
with  XBT  measurements  collected  by  the  R/V  LAMB.  DuTing 
the  acquisition  period  the  surface  conditions  were  moderate 
to  strong  - no  current  measurements  were  made.  The  data 
used  for  the  comparison  study  were  the  temperature  values  - 
uncorrected  and  corrected  for  vertical  displacements. 

Fig.  9 and  Table  7 give  the  results.  The  percent  dif- 
ference, %A  , between  the  XBT  and  buoy  data  are  given  in 
columns  4 and  6 - Table  7 - for  the  uncorrected  and 
corrected  measurements,  respectively.  The  VDCM  reduces 
the  by  a factor  of  three  for  the  given  conditions  of 

external  forcing. 

In  the  case  study  presented  above  the  XBT  measure- 
ment was  used  as  the  reference  assuming  no  instrumental 
error.  As  an  aside  the  BEAR  data  - corrected  for  vertical 
displacement  - were  compared  with  XBT  data  corrected  for 
i an  assumed  bias.  In  a study  conducted  for  the  U.S.  Coast 

Guard,  Gaillard  et  al  (1970)  showed  the  standard  sippican 
XBT  to  have  a response  bias  of  4.57  m or  2%  assuming  a drop 
rate  of  5.08  ms~^.  The  XBT  measurement  used  in  this  study 
was  corrected  for  the  above  bias  (given  in  column  7, 

Table  8)  and  compared  with  the  buoy  data.  The  reduction 
in  percent  difference  (column  8)  is  very  slight  - the 
order  of  1%  maximum. 


BEAR  12  Feb  76  1654  LST 


UNCORRECTEO  • 

CORRECTED  

LAMB  12  Feb  76  I403LST 

XBT  A 


Fig.  9 

Model  Results  vs  XBT  Comparison 
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Comparison  of  XBT  Data  - R/V  LAMB  12  Feb.  76 
1403  LST  and  BEAR  Buoy  1654  LST  - Uncorrected 
and  Corrected  for  Vertical  Displacement 
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4. 0 Sununary 

This  study  has  presented  the  mathematical  develop- 
ment of  the  mooring  ccmf iguration  model  (MCM)  and  the 
vertical  displacement  corrector  model  (VDCM) . 

The  MCM  uses  explicitly  the  surface  forces  (due  to 
wind,  wave,  and  current)  and  subsurface  drag  (due  to 
current)  to  compute  the  configuration  of  the  mooring. 

The  VDCM,  on  the  other  hand,  uses  the  total  measured 
surface  tension  and  approximates  the  subsurface  drag 
using  as  inputs  the  mooring  instrumentation  measure- 
ments - tension  and  line  inclination  at  the  surface  and 
inclination  at  the  cable,  nylon  junction.  The  MCM  is 
used  exclusively  to  analyze  mooring  configuration  for 
specific  case  studies  of  external  forcing.  The  VDCM 
is  used  to  correct  the  measured  temperature  values  for 
bias  due  to  vertical  displacements. 

Comparison  of  model  results  (MCM  vs  VDCM)  have 
shown  the  method  used  to  approximate  the  subsurface 
drag  to  be  valid  and  to  yield  results  consistent  with 
the  MCM.  Comparison  of  VDCM  vs  XBT  measurements  have 
shown  that  the  VDCM  reduces  the  error  between  buoy  and 
XBT  data  by  as  large  as  a factor  of  three  for  conditions 
of  moderate  to  strong  external  forces. 
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